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FINAL STAGES OF TRANSITION TO TURBULENCE IN PLANE CHANNEL FLOW 


By 

Sedat Biringen 

NASA/Langley Research Center (M.S. 359) 

Hampton, Virginia 23665 

This work, performed under Grant No. NAG-1 -228 from NASA/Langley Re- 
search Center, involves a numerical simulation of the final stages of 
transition to turbulence in plane channel flow. Three-dimensional, incom- 
pressible Navier-Stokes equations are numerically integrated to obtain the 
time-evolution of two- and three-dimensional finite-ampl itude disturbances. 
Computations are performed on the CYBER-203 vector processor for a 32x51x32 
grid. Results are presented for no-slip boundary conditions at the solid 
walls as well as for periodic suction-blowing to simulate active control of 
transition by mass transfer. Solutions indicate that the method is capable 
of simulating the complex character of vorticity dynamics during the various 
stages of transition and final breakdown. In particular, evidence points to 
the formation of a A-shape vortex and the subsequent system of horseshoe 
vortices inclined to the main flow direction as the main elements of trans- 
ition. Calculations involving periodic suction-blowing indicate that 
interference with a wave of suitable phase and amplitude reduces the distur- 
bance growth rates. 


1. INTRODUCTION 

Recent experiments by Nlshioka, Asia & II d a (1981) have shown that 
transition to turbulence In a plane channel flow follows a sequence of 
events similar to that observed by Klebanoff, Tidstrom & Sargent (1962) In 
the bound ary -layer transition. In this wrk, a direct nunerical integration 
of the Navier-Stokes equations is performed in an attempt to simulate these 
events in plane channel flow, during the later stages of transition. 

In their experiments, Nishioka et al . (1981) measured the streamwise 

i 

mean and fluctuating velocities, and ui, respectively, at a fixed 
streanwise location at a subcritical (linearly stable) Reynolds nimber, Re = 
5000 and simulated the various stages of transition by varying the distur- 
bance amplitude. Their observations show that subcritical instability takes 
place at a threshold ampl itude of (^i ) m -, v /^o = ^*01» where Ug is the mean 
velocity at the channel centerline. The evolution of this instability is 
evidenced by the intensification of the spanwise variation of the wavefront 
which develops into a peak -valley structure. Nishioka et al. (1981) 
observed that flow development follows a trend which is ‘similar to trans- 
ition in the boundary-1 ayer (Klebbnoff et al. 1962, Kovasznay, Komoda & 
Vasudeva 1962): local shear layers are formed away from the wall at span- 

i 

wise peak positions ( u x /U 0 = 0.06), these shear layers start to exhibit 
a "kink" which is the manifestation of secondary instability and is accom- 
panied by a "spike" in the u^Ug - 0.11). In rapid succession, two-, 
three-, five-, and multi-spike stages are observed with increasing ampli- 
tude of the primary disturbance. Nishioka et al. (1981) presented evidence 
that in the final stages of transition, the flow starts to develop struc- 


i. i 
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tures very similar to those found in fully developed wall turbulence. 

During this stage, the flow field is characteri zed by the development of a 
viscous sublayer, occurrence of the typical "streaks" close to the wall and 
the formation of horseshoe vortices sometimes referred to as the building 
blocks of wall turbulence (Theodorsen 1954). The present work simulates 
this sequence of events'. 

Effects of three dimensionality on transition have first been document- 
ed in detail by Klebanoff et al . (1962). Accordingly, three-dimensional ity 
manifests itself mainly in the spanwise velocity variations resulting in the 
production of streamwise vorticity which, in turn, interacts with the span- 
wise vorticity and drives the flow to breakdown. OrsMg & Kells (1980) and 
Patera & Orszag (1981) have expanded on this idea and studied the suscepti- 
bility of plans channel f-ow to three-dimensional vtu'Uirbances by nurierical- 
ly integrating the three-d rmensional Navier-Stokas equations. Their compu- 
tations at subcritical Reynolds nunbers revealed some interesting aspects of 
subcritical transition. They found that initially two-dimensional distur- 
bances which are finite-amplitude two-dimensional Orr-Sommerfeld eigen- 
solutions, decay slowly whereas the interaction of two-dimensional and 
three-dimensional disturbances drives the flow to instability for Reynolds 
ntmbers as low as 1000. In the present work we focus on the three-dimen- 
sional and nonlinear mechanisms that characterize transition in the final 
stages and during flow breakdown. Hence, we employ the full three-dimen- 
sional time-dependent Navier-Stokes equations. Calculations were performed 
at a linearly stable Reynolds nunber (Re = 1500), with finite-ampl itude two- 
and three-dimensional eigensolutions of the Orr-Somerfeld equation used as 
the initial conditions. No attempt was made^herein to study the effect of 


different Initial conditions or of Reynolds nunbers; this Is the subject of 
another investigation. We also Investigated the Idea of using periodic suc- 
tibn-blowing at the solid walls as a means of controlling transition. These 
computations were performed at Re * 7500 which Is the linearly unstable 
range for this flow. Since transition control must be applied prior to the 
flow breakdown (e.g. before the occurrence of multi-spikes), disturbance 
amplitudes during these calculations were not large enough to Impose stiff 
restrictions on the allowable time step. 

In section 2, of this report, the nunerical methods used in this study 
are briefly discussed. In section 3, results of calculations are presented 
and finally in section 4 a summary of results and some concluding remarks 
are given. 


2. THE CALCULATION PROCEDURE 


The calculation procedure Is based on the Incompressible Navier-Stokes 
equations In primitive-variable form, 


3u 


_!+«.!!!- -1 lEl + v 

at * ax^ p ax 1 a \ 3x ^ 


and the continuity equation, 
au 4 

L - n 


0) 

( 2 ) 


where u^ are the velocities along the x^ directions, p Is the density, 
v is the kinematic viscosity and p' is the total hydrostatic pressure. 

The equations are non-dimensional 1 zed by the mean centerline velocity Uo 
and the channel half-width, h. The flow is assuned to be driven by a con- 


stant mean pressure gradient 2/Re, where Re is the Reynolds number 
given by U 0 h/v. Also, the convective terms are written in a form which 
prevents occurrence of nonlinear instability in the numerical solution 
procedure by ensuring conservation of momentum and energy (Mansour, Per- 

i 

ziger, & Reynolds 1978). The final form of the Navier-Stokes equations 
reads, 


3 U . /3U. 3U„ 
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where P = p'/p + i u^u^ is the pressure head and 6^ 


is the Kronecker 


delta. 


The flow 1c assuned to be periodic In the streamwlse and the span- 
wise x 3 directions along which the flow field variables can be expanded In 
terms of Fourier series. This enables the use of the pseudo- spectral method 
(Orszag 1972) to calculate the spatial derivatives along xi and X 3 by 
use of discrete Fourier transforms. Considering transforms in the x^ direc- 
tion, along which there are Nj. equally spaced mesh points, the velocity 
component u^ can be written as 

N/2 -1 ^ y 

ui(xi) * z ui(ki)e 1kl 1 (4) 


where Xi B maxi , m = 0,1,..., N-l and ki = 2-rrni /Wi axi . Accordingly, the 
Fourier transform of u x is 

= i V m(x,).- ,k * x > (5) 

Ni m=0 


The spatial derivative of 

Ui along xi 

can now be written as 

3Ui (xi) 

Ni/2-1 
= s 

ikiui (kj )e lklXl 

3Xi 

ni= -Ni/2 



The derivative can be computed by forming the Fourier transform of m (xi), 
multiplying the result by ik 1 and computing the inverse transform. For 


periodic functions, the pseudo-spectral method provides a means by wh'ch the 
spatial derivatives are evaluated with maximun accuracy for a given nimber 
of grid points. Along the x 2 -d1rect1on, a mesh stretching that concen- 
trates grid points close to the solid walls Is employed (Moln, Reynolds & 
Firziger 1973). The resulting mesh enables the resolution of the sublayer 
that is formed during transition fcr y+ < 2 , where y+ Is the coordinate 
along x 2 , Inwall units. Spatial derivatives along x 2 are evaluated by 
a second-order finite-difference scheme on this non-uniform mesh. 

The governing equations were mmerlcally Integrated by the semi-implic- 
it method of Moln et al. (1978). This procedure employs the explicit Adams- 
Bashforth method for the convective terms and the implicit Crank -Nicholson 
method for pressure and for the viscous diffusion terms. In order to start 
the two time-level Adams-Bashforth method, the Euler- implicit method is used 
at the first time step. 

Once the governing equations are discretized in time, a two-dimensional 
Fourier transform along the periodic directions xi and X 3 transforms the 
equations into the k 1 -k 3 wave-nunber space. The transformed equations are 
written below in block-tridiagonal form for inversion along X 2 


A F^J + B f!J +1 

— ~j+l — j 


+ C F n , + ] 
- ~J-1 


R ? 


( 7 ) 


In (7), A, B, and C_ are coefficient matrices, Fj is the solution vector 

at the advanced time level, n+ 1 , and at the X 2 -directional node, j; fij is 
right-hand side vector that contains the convective, diffusive and pressure 
terms at the previous time levels. These are given as 
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Coefficients of the. finite-difference operators that appear In the matrices 
A, B, and C are given as 


fi * ^j + i<V + V> B 2 j = - 2 /A j-i A j* A 2 J = 2 / VVi + V 

Clj - -Alj = l/(A j+1 - Aj) 

I 

A j+1 = ^ X2 ^j+1 “ ( X2 )j 

Since all the flow variables in the solution vector contain an imaginary and 

* 

a real part, the block-inversion process is applied twice for each pair of 

kj and k3, which are the waventmbers along xi and X3, respectively. 

The assumption of periodicity in x* and x 3 eliminates the^necessity 
of applying explicit boundary conditions along these directions. However, 
due to the presence of solid boundaries along the x 2 direction, no-slip 
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boundary conditions are imposed on u x , u 2 , and u 3 and the pressure at 
he wall Is calculated by a second order approximation from the interior of 
;he flow field. That the pressure boundary conditions are consistent with 
the x 2 -momentum equation at the wall, has been shown by Moin et al. (1978). 

Initial conditions were prescribed from the two- and three-dimensional 
eigensolutions of the Orr-Sommerfeld equation by considering that even for 
subcritical Reynolds numbers, plane channel flow can be driven to instabili- 
ty if tho least stable two-dimensional finite-amplitude Orr-Sommerfeld 
eigenmodes are interacted with f inite-ampl itude three-dimensional eigenmodes 
(Orszag & Kells 1981). The most explosive situation arises when the three- 
dimensional eigenmodes are aligned with the main flow direction at ±45 to 
±60 degrees. Accordingly, we have used the following initial condition: 

u(x) * U(x ? ,0,0) + u 2D (x 2 )e laX i + u 30 (x 2 )e laX i ±m 3 (8) 

Here, U(x 2 ,0,0) is the parabolic velocity profile of plane channel flow. 
The eigenfunctions u 20 (x 2 ) anc * U 3D^ 2 ^ correspond to two-dimensional and 
three-dimensional solutions of the Orr-Sommerfeld equation at Re = 1500, 
respectively. The two-dimensional solution was obtained for a - 1 whereas, 
the three-dimensional solution was obtained for a - 1, p - ±1. A computer 
program, which essentially uses the Kaplan filtering technique was used for 
the solution of the Orr-Sommerfeld equation (Reynolds 1967). The final 
amplitudes were chosen so that the maximum value of the Xj-directional two- 
dimensional disturbance was set equal to 0.1 1U 0 and the maximum amplitudes 
of the x^directional three-dimensional disturbances were each set equal to 
0.05U 0 . 


3. RESULTS AND DISCUSSION 


The finite-difference system (7) was solved on the CYBER-203 vector 
processor at NASA/Langley Research Center. A 32x51x32 mesh was employed 
along the y>i~, X 2 -, and X 3 -directions, respectively. The computer code was 
fully vectorized and vectorized library subroutines were used for the main 
computational operations that the solution technique anploys. These vector 
operations are mainly one-dimensional fast Fourier transform (FFT) to cal- 
culate spatial derivatives with the pseudo-spectral method, two-dimensional 
FFT to transform the equations into k 1 -k 3 wave-number space and block- 

tridiagonal matrix inversion along X 2 . For the FFT operations, typical 
vector lengths were around 1000 , which is an optimal vector length to take 
full advantage of the vector processor. For the block-trid agonal matrix 
inversion (which essentially is a scalar operation), a vectorized subroutine 
that inverts a large nunber of tridiagonal systems simultaneously was used. 

This procedure decreases CPU time significantly by reducing the number of 
scalar operations required to invert each system separately. Tne fully 
vectorized code takes about 30 sec. of CPU time per time step for the 32x51x32 
mesh to solve the finite-difference system (7) on a computational box, in which 
the flow is confined between rigid walls at x 2 = ±1. Periodicity lengths 
(box 1 engths) along xi and X 3 were chosen so that the smallest wave 
nimbers allowed in the computational domain were equal to a = 1 and 3 = 1 , 
respectively, i.e., the box length was set equal to 2 ir along these directions. 

It should be recalled that the time- advancement scheme employed in this 
work is partly explicit (on the convective terms) and partly implicit (on 
the diffusion and pressure terms). Although. in view of linear stability 
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analysis, implicit methods are unconditionally stable (extrapolation to 
nonlinear equations is sometime vague), the mixed nature of the present 
scheme as well as the time-accurate nature of the problem under investiga- 
tion necessitate adherence to stability bounds of explicit schemes. There- 
fore, in all the calculations reported herein, the convective stability 
condition (the Courant-Friedrichs-Lewy condition) that requires the Courant 
nimber (C.N.) to be always less than one and the diffusive stability condi- 
tion ware obeyed. With (Ax 2 ) mi - n = 0.0092 and at = 0.025, where Aj j S the 
nondimensional time-step, through the course of the calculations C.N. varied 
as 


C.N. = AT 


ui 

AXl 


U2_ 

AX2 


U3 

AX3 


< 0.2 


( 9 ) 


rmax 


whereas the diffusive stability criterion, D, varied as 

AT 


D = L 


Re 


(A x) ' 


mm J 


< 0.04 


( 10 ) 


so that the diffusive stability condition which requires D <0.5, was also 
always satisfied. The computer program was tested by calculating the growth 
rates of smal 1-ampl itude Orr-Sommerfeld waves. For a wide range of Reynolds 
numbers (between 1000 and 10,000) the agreement between the computed re- 
sults and the linear theory was better than 0.5%. 

In the subsequent parts of this section, results obtained from the 
numerical integration of the finite-difference system (7) for the time-evo- 
lution of the initial disturbances are compared with the experiments of 
Nishicka et al . (1981). It should be noted that this experiment was done at 
Re = 5000 (subcritical) whereas the computation was done at a lower 
aubcritical Reynolds number. Re = 1500. The selection of a higher Reynolds 
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nunber {e.g., in the linearly unstable range) makes the governing equations 
very stiff and requires the use of extremely small time steps for nunerical 
stability. Hence, the selection of Re = 1500 was mainly to force the 
computations Into transition and breakdown with the least amount of computer 
expense. It should be noted that wall phenomena characteristic of the final 
stages of transition are essentially independent of Reynolds nunber as 
observed by Nishioka et al. (1981) for channel flow and by Smith & Metzler 
(1982) for boundary-layer transition. Therefore, the difference between the 
Reynolds numbers of the experiment and the simulation should not have any 
important consequences for the qualitative comparisons between the two sets 
of results. 

3a. Plane Channel Flow with No Mass Transfer 

In their experiments, Nishioka et al . (1981) identified the various 
stages of transition according to the number of spikes appearing. in the os- 

t 

cilloscope traces of the disturbance velocity, u^ Since the time-axis of 
the experiment is interchangeable with the xi-axis of the computation, we 

i 

obtain similar traces by plotting u x along Xj (over two periods) as shown 
in figure 1. In this figure, the first frame shows the sinusoidal variation 
of the initial conditions at T = 0.. This is followed by the nonlineaer dis- 

t* ' 

I 

tortions of the initial conditions resulting in variations of Ui which 
strongly resemble the oscilloscope traces at the one-, three-, and five- 

I 

spike stages of the laboratory flow. Ii. particular, the variations of ui 
with x L at T = 44 are very similar to the ensemble- averaged waveforms 
presented by Nishioka et al. (1981) at the five-spike stage. The last frame 
in figure 1 shows variations of at a time (T - 79) much later than the 

occurrence of the five-spike stage (T = 44). Hence, we note the absence of 


flotations of rich frequency content that are characteristic of fully- 
developed turbulence; this lends support to the Idea that the last stage of 
transition may not be spontaneous and explosive a phenomenon as generally 

supposed. However, at this stage, the u| variations are very similar to 
velocity oscillations observed in wave packets and indicate the occurrence 
of "patches" of turbulent fluid. For the remainder of this discussion, 
results from the simulation at the spike stages (shown in figure 1) will be 
used, whenever possible, for comparisons with the corresponding spike stages 
of the experiments of Nishioka et al. (1981). 

In figure 2, we show a history of the time-evolution of the flow in 
terms of the maximum amplitude of the two dimensional primary disturbance, 
its two-dimensional harmonic and the three-dimensional primary disturbance. 
The trends displayed by these quantities are generally similar to the re- 
sults of Orszag & Kells (1980) which they obtained from computations per- 
formed at Re = 1250. The main features of these trends are the rapid de- 
crease in the two-dimensional primary-wave amplitude and the rapid increase 
of the amplitude of its harmonic. ’Also, the three-dimensional primary-wave 
amplitude first increases, then gradually decreases at around T = 30. In 
the present calculations, we observe that variations of the amplitudes start 
to fluctuate as early as T * 15 but, even at later times, the fluctua- 
tions do not display an explosive trend. That no "explosive" instabilities 
were found in the present computations is in accord with the findings of 
Nishioka, Asai & Iida (1980), which imply that breakdown in channel flow is 
as gradual as the growth of instabilities found in free shear flows. 

In figure 3, we present plots of maximum root-mean-square (rms) ampli- 


tudes of ui, (ui over two periods along X 3 . Here, we define the 

rms value as an average over Xj, . We note that at T - 0, the wave pat- 
tern is sinusoidal and corresponds to peaks at X 3 a 0 , 2 ", W and to 
valleys at x 3 - u, 3rr. Subsequently, at later times the nonlinear dis- 
tortions of the wave front result in minima occurring at the peaks and 
maxima occurring at the valleys in accordance with the experiments of 
Nishioka et al. (1980). At later times, an Increase in the frequency of the 
peak-valley structures is clearly evident suggesting an increase in the nun- 

ber of characteristic vortex structures along X3. 

Plots of velocity profiles averaged over the X!-x 3 plane, <u 1 >, are 
shown in figure 4 for laminar (initial), late transition and "early turbu- 
lence" stages at T = 0, T = 44 and T = 79, respectively . The <u^> 
distribution at T = 79 has a strong resemblance to the turbulent 
channel flow profile, with increased velocity gradient at the wall and with 

a full profile indicative of turbulent mixing. Although, as expected, <ui> 
profiles do not show any fluctuations (or inflections), plots of instantan- 
eous velocity profiles do show very* strong inflections, especially in the 
regions close to the walls (figure 5). This indicates that the interac- 
tion of two- and three-dimensional waves close to the walls is the central 
mechanism that drives the flow to instability. This is in accord with the 
idea that the flow will undergo transition only for a selected band of span- 
wise wave numbers, the most "dangerous" of which result in three-dimensional 
disturbances with maxima occurring close to the walls (Orszag & Kells 1980). 
In figure 6 , the velocity profile <u+^ = <ui>/u T versus y+ = X 2 i^/v is 
plotted; here u T is the friction velocity and is calculated from d<ui>/dx 2 


at the wall times 1/Re. Although, the plots indicate the formation of a 
sublayer, and the change from T *44 to T * 79 shows a gradual approach 
to the law-of-the-wal 1, the difference is still especially apparent in the 

I 

logarithmic region. At T = 44, the Reynolds nimber based on friction 
velocity is equal to 69 and, as expected, is larger than its Initial 
(laminar) value of 55. 

Plots of pi ane- averaged fluctuations intensities, <{ui~ <Ui>)2>, 
are shown in figure 7 at various T. There are several interesting features 
of this figure. Firstly, at T - 17, the shift in the position of peak 
amplitude towards the channel center (X 2 K -0.6), as well as the increase 
in the maximum amplitude indicate that the development of the computed flow 
field is compatible with experimental observations pertaining to the one- 
spike stage of the transition process (Tani 1969). It will later be shown 
that during this stage there is a substantial increase in the spanwise 
vorticity, a> z , away from the lower wall around x 2 = - 0.6. Secondly, in 
accordance with the laboratory flow of Nishioka et al. (1981) at later 
stages, the computed intensity profile displays a second peak occurring 
close to the wall associated with turbulence production. At T = 44 cor- 
responding to the five-spike stage, we see that peak intensity has reached a 
value typical of turbulent channel flow; however, the peak occurs uncharac- 
teristically away from the wall. Finally at T = 79, the peak in the inten- 
sity profile has moved towards the wall but even at this stage, the distri- 
bution has not assuned the asymptotic f ully-turbulent form. In figure 8, 
plots of pi ane- averaged shear stress, <ui - <ui>)u 2 >, profiles are shown. 
The increase in the magnitude of shear stress from T - 0 to T = 17 clearly 
indicates the effects of nonlinearity in transferring energy from the mean 


flow to the fluctuating motion. At the five-spike stage (T a 44), the 
maximum shear stress has attained a value typical of turbulent channel flow. 
However, the locatiion of the maximum is away from the wall and roughly cor- 
responds to the location of the peak in the respective intensity profile 
(figure 7), It can thus be inferred that at this stage of the comutation, 
the energy-exchange mechanisms of fully-developed wall turbulence is not yet 
reflected by the p 1 ane- averaged velocity correlations. 

Spanwise variations of Uj are plotted at various distances along x 2 
in figures 9a-9e. In figure 9a, we show initial distributions which are 
highly three-dimensional. At T = 17 (figure 9b), the distributions indicate 
a velocity defect both at peak and valley as well as a velocity excess in 
between. This is indicative of intensification of the initial streamwise 


vortices and appearance of weaker streamwise vortex pairs in accordance with 
the Bermey-Lin (1960) theory. Subsequently, figures 9c-9e display the form- 
ation of additional vortex pairs indicating transport of energy down the 
wave-nimher spectrum. The spanwise synmetry imposed by the initial condi- 
tions is retained through the five-spike stage. At this stage, an estimate 
of the flow field resolution can be obtained from the spanwise distance be- 
tween peak positions, \ z . Nondimensional ised by u 2 and v f typically * z = 
115. This is larger than but comparable to X = 80, which is the typical 
spanwise length in the laboratory flow during the five- spike stage (Nishioka 
et al. 1981). It should be noted that the spanwise characteristic length in 
wall turbulence is about 100. Therefore, it could be asserted that at this 
stage present results are representative of initial wall turbulence. 

A more detailed description of the transition process can be obtained 
from contour plots of equi-shear lines, 3 Ui/ 3X2 (which correspond to 
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approximate spanwise vorticity, w z ) in the x x - x 2 plane at the position of 
maximum u^Uq, Results from the computation that correspond to the various 
stages of the laboratory flow are presented in figures 10-15 between the 
lower wall (x 2 15 -1.0) and channel center (x 2 - 0.0). In figures 11a - 
13a, figures 4-6 of Nishioka et al. (1981) are also shown for a qualitative 
comparison with the present results. In figure 10, contour plots corre- 
sponding to the initial conditions and In figure 11, contour plots corre- 
sponding to the "one spike'* stage (T 3 0 and T 3 17, respectively) are 
shown. In both the laboratory flow and the computation, the typical head of 
the shear layer appears very clearly at the one-spike stage, Indicating the 
formation of a shear layer away from the wall at about x 2 3 - 0.6 due to the 
induced velocity from the streamwise vortex system. In addition, the sud- 
den dip of the shear layer from the high-velocity outer flow to the low- 
velocity region clearly appears as a kink in both the computation and the 
laboratory flow. Since the grid points are finely clustered along x 2 
close to the wall, vorticity concentrations in this region are also ade- 
quately resolved by the numerical simulation. 

Figure 12 shows equi -shear lines at T = 27 corresponding to the three- 
spike stage of the laboratory flow. Due to the secondary instability mani- 
fested in the previous stage, breakdown of flow structures into smaller 
scales is observed in both the experiment and the computation. The growth 
of the kinked portions of the equi-shear lines into the so-called "hairpin 
eddies" is clearly depicted in the experiment. The computation displays a 
similar evolution: the head of the shear layer is lifted up towards the 
channel centerline, and the kink in the shear layer is quite apparent in 
this stage of the simulation. Simultaneously witth this activity taking 


:lace in the outer (high speed) portions of the flow field, both the 
experiment and the computation show an Intense shear layer developing close 
to the wall, which Is Indicative of turbulence generation. It Is generally 
agreed that hairpin eddies which are lifted towards the centerline erupt 
into turbulent spots. However, figure 12 indicates that wall turbulence 
may also be closely associated with vorticity dynamics simultaneously taking 
place with the eruption of hairpin eddies. Contours of equl-shear lines at 
T - 44 corresponding to the five-spike stage of the experiment are shown in 
figure 13. In both the experiment and the computation, the intense shear 
layer developed in the wall region is closely discernible. It should also 
be noted that, the spanwise position of maximum w z which appears very 
close to the wall and does not necessarily coincide with the spanwise 
location of maximum u^Ug where the contours are presented. The most sig- 
nificant feature of figure 13 is the existence of distinct vortex structures 
in the wall region both in the laboratory flow and in the computation. 

These vortices are inclined to the main flow direction at an angle that 
varies between approximately 14° to 40° and show a close resemblence to the 
energetic horseshoe vortices characteristic of initial wall turbulence (Hama 
& Nutant 1963, Klebanoff et al, 1962). It is these vortices roughly aligned 
along the direction of maximum extensional stress that are mainly responsi- 
ble for extracting energy from the mean shear (Tennekes & Lumley 1972). 

In figure 14, equi-shear lines at spanwise locations x 3 = (x 3 ) 0 + 2ax 3 
and x 3 = (x 3 ) 0 + 4Ax 3 at T « 44 are compared. Here, we define (x 3 ) 0 as 
the spanwise position corresponding to figure 13 and Ax 3 is the mesh size 
along x 3 . Figure 14 suggests that the shear layer is formed from a system 
of horseshoe vortices in succession such that as the first-born vortex 
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erupts Into the channel center, a new one forms at the wall. As evidenced 
In Figure 5, the vortex lift-up Is still clearly discernible at the "early 
turbulence" stage (T = 79). Here, we note the vortex lift-up towards the 
channel center but Its extension Is considerably shorter than that observed 
previously at the five-spike stage. This variance suggests the existence of 
different vortex structures characteristic of these two stages. In fact, we 
shall later show that the typical horseshoe vortex of 1 ;te transition Is not 
the basic structure of early turbulence. 

As we have noted before, the mechanism that is responsible for the 
generation of vortlcity concentrations Is usually explained as vortex 
stretching (and deformation) by the mean flow. The stretching and deformed 
layer moves downstream with a translation velocity that induces lower local 
velocity of the upstream edge of the vorticity layer than Its downstream 
edge (Komoda 1967). An examination of m contours of (figures 16 and 17) 
along with the approximate spanwlse vortlcity contours at T = 27 (figure 
12b) and at T = 44 (figure 13b) in the x* - xg plane clearly shows a similar 
trend. We observe that the nose of the vorticlty layer is general ly associ- 
ated with higher velocities than the upstream region and large variations of 
the local velocity exist within the layer. 

Finally, normal velocity, u 2 , contours are shown in figures 18-21 
during the three-spike, multi-spike and early-turbulence stages at the 
position of maximum u x /U 0 . Figure 26 shows u 2 -contours in the X! - Xa 

plane at T = 27, correspond! ng to the three-spike stage of the laboratory 
flow. In this figure, there is clear evidence of the beginning of an alter- 
nating up-and-down flow similar to the pattern described by Kovaszny et al. 
(1967), which is characterized by the intense updrift accompanied by fluid 


drifting down at both sides. This alternating structure of up-and-down flow 

in the Xi - x 2 plane is more evident at T B 44 (figure 19), corresponding to 
the multi-spike stage of the experiment. Also In this figure, alternating 
regions of fluid with scales smaller than those at T * 27 are depicted. In 
figure 20 and 21, normal velocity contours are shown In the x 2 - x 3 plane at 
T « 44 and T = 79 corresponding to the multi-spike and early-turbulence 
stages, respectively. These contours are plotted in the region between the 
lower wall (x 2 3 -1,0) and x 2 n -0.92, display the evolution of alternating 
structures very similar to the characteristic streak-like structures found 
In the wall region of turbulent channel flow (Moin & Kin 1982). 

3b. Active Control of Transition by Periodic Suction-Blowing 

In this section we demonstrate the usefulness uf employing periodic 
suet ion- blowing boundary conditions at the solid walls for controlling 
transition. All these calculations were performed at a Reynolds nimber well 
within the linearly unstable range, i.e. at Re a 7500. The main idea here 
Is to cancel or modify the “most dangerous" disturbance by interfering with 
a wave of suitable amplitude and phase. This idea has been proven to be 
effective in controlling boundary-layer transition in an experiment in which 
the two-dimensional Tollmien-Schlichting wave was effectively cancelled by a 
control wave generated downstream of the initial disturbance. The amplitude 
and phase of the control wave (produced by a second wave generator placed 
downstream of the first) was adjusted to minimize the maximum disturbance 
from the first wave generator (Milling 1981). Expanding on this idea and 
assuming the two-dimensional primary wave to be the most dangerous disturb- 
ance, in this work we prescribe the boundary conditions at the solid walls 
such that 


(a) the resulting control wave is compatible with flow periodicity 
along x 1 and x 3 so that there will be no net mass transfer into 
the flow field; 

(b) the resulting control wave will be 180’ out of phase with the two- 
dimensional primary wave; 

(c) the control wave amplitude is adjusted for maximum attenuation of 
the primary disturbance. 

In figure 22, results for a case in which the Initial disturbances are 
two-dimensional solutions of the Orr-Sommerfeld equation are shown. The 
maximum amplitude of u L was set equal to 1 % of the channel centerline 
velocity. Periodic suction-blowing boundary conditions were prescribed at 
each wall at T - 5.3 for one time stop, before and after which the channel 
walls were made impermeable via the no-slip boundary conditions. The most 
effective wave produced by the suction-blowing boundary conditions was found 

i 

to be 180’ out of phase with the two-dimensional primary u 2 (corresponding 
to O 2 (l,0) in the Fourier space) with twice its amplitude. The effect of 

i 

the control wave to reduce the growth rates of both (u x ) ^-primary and its 
harmonic is evident from figure 22. In fact, the control wave has resulted 
in a decaying harmonic preventing energy transfer down the wave-nunber 
spectrum, consequently inhibiting the. proper development of the transition 
process. 

To test the effects of three-dimensionality a similar calculation was 

* 

performed in which three-dimensional solutions of the Orr-Sommerfeld equa- 
tion were also included into the initial disturbance field. For this case, 
the maximun amplitudes of the two- and three-dimensional disturbances were 


i « 
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set at 7 >% and of the channel centerline velocity, respectively. The most 
effective control wave was found to be 180* out of phase with ( u 2 ) g D with 
2,5 times Its amplitude. The characteristics of the resulting flow field 
are shown in figure 23, after the periodic suction-blowing boundary 
conditions are applied at T * 5.3 for one time step. In sunmary, the 

i i 

control wave attenuated the growth rates of (ui) total and ( u i ) g D “ primary, 

* 

and resulted in a decaying (*Ji ) go " harmon i c as in the previous case. Note 

i 

that ( u i ) 3 d~ Primary remained unchanged. 

It is apparent from the preceding discussion that the application of 
periodic suction-blowing boundary conditions on the solid walls of the 
channel has a considerable effect in reducing the growth rates of two- and 
three-dimensional disturbance fields. In practice the technique can be im- 
plemented via a suction-blowing strip along the body span. The spectrun of 
the disturbance field just upstream of the strip must be available in order 
to activate the correct boundary conditions and obtain the most effective 
control wave. Further work is necessary to evaluate the frequency-response 
and feedback delay of such a mechanism. 
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4. SUMMARY AND CONCLUDING REMARKS 

In this study, final stages of transition to turbulence In plane 
channel flow have been simulated by a direct nunerlcal solution of the 
Navier-Stokes equation It is found that, In spite of the limited resolu- 
tion of the 32x51x32 grid employed in the computations, the simulation is 
capable of reproducing most of the essential features of wall phenomena 

observed In the laboratory. Grid resolution in the xi and X 3 directions, 

along which the flow is periodic, is found to be adequate to capture the 

sequence of events that lead to early turbulence. Vorticity contours in the 

vicinity of the lower wall indicate formation of a system of horseshoe 

vortices with legs or extensions in the xj, - X 3 plane comprised of 
counter-rotating streamwise vortex pairs. Our findings are also in accord 
with the Benney-Lin (i960) theory, e.g., the present computations clearly 
depict frequency-doubling of streamwise vorticity. It should be noted that 
the initial conditions used in this work are of the Benney-Lin type and 
consist of a vorticity field with a strong spanwise component and weak 
streamwise and transverse components. Therefore, the question of the origin ’ 
of stream-wise vorticity and its precise relation to distortions of spanwise 
vorticity remain to be addressed in future work. 

At later stages of the computation, transverse velocity contours indi- 
cate the formation of streak-like structures alternating in the spanwise 
direction. Typically, the spanwise characteristic length of these vortices, 
inferred from the spanwise variations of u^ was found to be approxi amtely 
X 7 = 115; this is close to X 2 - 100 of fully-developed wall turbulence. 

It was found that during the later stages of transition, flow field 
statistics indicate the formation of a laminar sublayer; however, the 


development of the logarithmic region and consequently the approach to 
fully-developed turbulence Is slow. This gradual approach to steady state 
Is also reflected in the profiles of pi ane- averaged Intensity and shear 
stress. 

The mein deficiency of this study stems Inevitably from limited spatial 
resolution and manifests Itself in several ways. Firstly, at later stages 
of the computation, insufficient mesh resolution results in lower gradients 
of the mean velocity in the viscous sublayer. This, in turn, causes less 
turbulence production in the wall region as evidenced by the absence of 
peaks in the intensity and shear-stress profiles close to the wall. Second- 
ly, the finite cut-off wave nunbers along x x and x 3 preevent the formation 
of a proper wave-nimber spectrum. As a consequence, pollution of Fourier 
modes and accumul at ion of excess energy at low wave numbers become very 
significant sources of accuracy at large T. Hence, proper simulation of 
transition beyond the early-turbul ence stage necessitates the use of higher 
grid resolution. Even then, the incorporation of a mechanism to account for 
the subgrid scale turbulence would be required for a more accurate and 
realistic representation of the flow field phenomena. 

This work was supported by NASA/Langley Research Center under Grant 
No. NAG-1 -228. The author is indebted to P. J. Bobbitt and W. D. Harvey for 
their interest and encouragement during the course of this work. 
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2. Time-history of maximum disturbance amplitudes. 

(a) Two-dimensional primary; (b) three-dimensional 
primary; (c) two-dimensional harmonic. 
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Figure 4. Plots of planer-averaged mean velocity profiles 
(a) T«0; (b) T=44; (c) T= 79. 



Figure 5. Plots of instantaneous velocity profiles at x^ 
xo= tt/ 2. (a) T=17; (b) T=44; (c) T= 79. 








Figure 6 


Plots of plane-averaged mean velocity in wall units. 

A, T - 79. 



Figure 7. Plots of plane-averaged fluctuation intensities of u 
(b) T = 17; (c) T = 27; (d) T - 44; <e) T = 79. 
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Figure 12. Contour plots of 3u 1 /3x 2 . (a) Three-spike stage, figure 5 from 

Nishioka et al. (1981); (b) computations at T » 27 in the xi-x 2 
plane, contours from 0.3 to 6.8. 





Figure 13. Contour plots- of Su^/Sxn. (a) five-spike stage, figure 
6 from Nishioka et al. (l981) ; (b) computations at T^A 
in the Xj-X 2 plane; contours from 0.7 to 7.0. At this 
spanwise position, we define xoE(x 3 ) . 
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Figure 16. Contour plots of Uj in the x^-x 2 plane at T = 27. 
from 0.0 to 1.08. 
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Figure 22. Two-dimensional wave cancellation 





Figure 23. Three-dimensional wave cancellation 





